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La pr£sente invention a trait k une methode d' evaluation d'une courbe de pression 
capillaire des roches d'un gisement souterrain k partir de mesures sur des debris de rpche 
qui y sont preleves 

10 Etat de la technique ) 

Mesures en laboratoire sur carottes ou debris de forage .> 

La mesure de paramStres p£trophysiques tels que la permeabilite, la porosite et les 
propri6t6s capillaires sur des fragments de roche remontes au cours du forage d'un puits.Au 
travers d'un formation souterraine, constitue une opportunite interessante pour les 
15 compagnies op£ratrices d'obtenir rapidement une premi&re caract6risation pdtrophysique 
de zones productrices traversees par le puits. 

Par le brevet FR 2 809 821 du demandeur, on connait un systeme pour evaluer des 
parametres physiques tels que leur permeability absolue de roches poreuses d'une zone 
d'un gisement souterrain, k partir de d6bris rocheux remontes dans la boue d'un forage. 

20 Dans une enceinte ou les debris sont plong6s dans un fluide visqueux, on injecte de ce 
mfeme fluide sous une pression croissante avec le temps, jusqu'k un seuil de pression 
definie, de mani&re k comprimer le gaz pi6g6 dans les pores de la roche. Cette phase 
d'injection est suivie d'une phase de relaxation avec arret de l'injection. L* Evolution de la 
pression durant le processus d'injection ayant et6 mod61is6e a partir de valeurs initiales 

25 choisies pour les parametres physiques des debris, un calculateur les ajuste it6rativement 
pour faire coincider au mieux la courbe de pression mod61is6e avec la courbe de pression 
r6ellement mesurSe 



Par la demande de brevet FR 02/0023 du demandeur, on connait une autre methode 
pour evaluer des parametres physiques tels que la permeabilite absolue et la porosite des 
roches d'une zone d'un gisement souterrain, a partir e-galement de debris de forage. Une 
enceinte contenant les fragments de roche et remplie d'un fluide visqueux est nnse en 
5 communication avec un reservoir de ce meme fluide sous une pression definie de mamere 

a comprimer le gaz piege dans les pores de la roche. La duree duplicat ion cfecette ^ 
"~^Tn^n"qu'elle est courte ou plus longue permet de mesurer sou la vanafibli Arhr - 
pression dans 1'enceinte ou la variation du volume effectivement absorbe par les fragments 
de roche. Ensuite, on modelise Evolution de la pression ou du volume dans l'encemte, a 
0 partir de valeurs initiales choisies pour les parametres physiques des fragments, et on ajuste 
iterativement les valeurs des parametres physiques des fragments de roche pour que 
revolution modelisee s'ajuste au mieux avec Involution mesuree du parametre physxque 
dans F enceinte. 

Dans le domaine de la caracterisation petrophysique, la pression capillaire est aussi 
15 une donnee tres importante pour les operateurs car elle conditions : 

e la distribution initiate des fluides dans le reservoir depuis la zone aquifere (dite 
WOC pour "Water Oil Contact" par les gens de Tart) jusque dans la partie superieure du 
reservoir (zone de transition). Suivant la courbe de pression capillaire associ6e a une roche 
r6servoir et la nature des fluides en place, cette zone de transition pent s'etendre sur 
20 quelques metres ou quelques dizaines de metres ce qui a un impact important sur la 
determination des accumulations en place. 

s la pression d'entree d'une roche ce qui est parti culierement important pour les 
roches couvertures. Par exemple, pour un reservoir de stockage de gaz, la pression d'entree 
des roches de couverture conditions directement la surpression admissible dans les 
25 niveaux de stockage sans avoir de fuites. 

Avec les techniques actuelles, la courbe de pression capillaire est obtenue par 
Tintermediaire de mesures en laboratoire sur des carottes de gisement. Ces methodes sont 
cheres aussi bien a cause des operations le carottage que des operations de mesure sur les 
carottes et les resultats ne sont souvent disponibles que plusieurs mois apres le forage. 
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Approches pour obtenir rapidernent la courbe de pression capillaire 

II existe cependant des methodes alternatives decrites dans la litterature pour 
evaluer le courbe de pression capillaire rapidernent, soit en cours du forage soit en leger 
differe. 

5 L' approche la plus couramment utilisee consiste k utiliser la technique de la 

porosimetrie au mercure pour mesurer la courbe de pression capillaire Pc air/mercure 
directement & partir de d^blais de forage ou "cuttings". Neanmoins, la courbe obtenue 
differe largement par rapport a la courbe de reference obtenue sur carotte aux fortes 
saturations en fluide mouillant. D'autre part, cette approche repose sur Tutilisation de 
10 mercure qui est extrSmement polluant et progressivement interdit par la legislation dans de 
nombreux pays, ce qui pose un probleme majeur pour 1' application de cette technique dans 
un futur proche. 

Une autre methode connue utilise la Resonance Magnetique Nucleaire (RMN ) 
pour estimer rapidernent la courbe de pression capillaire a partir des donnfies diagraphiques 
15 mesurees dans le puits peu apres le forage. Elle est mise en oeuvre notamment dans les 
publications suivantes : 

- . Bowers, M., A. et al. : " Prediction of permeability from capillary pressure curves 

derived with NMR", 17 September^ 1 998 j 5 

- Marshall, D., et al. -."Method for correlating NMR relaxometry and mercury injection 
20 data" , SCA n° Society of Core Analysts Internationa] Symposium 1995 ; 

- Volokitin, Y., W. J. et al. :"A practical approach to obtain 1st drainage capillary 
pressure curves from NMR core and Log data" , SCA n° Society of Core Analysts 
International Symposium 1999. 

Le signal de relaxation RMN est d'abord converti en termes de distribution de taille 
25 de pore puis ensuite en terme de distribution de taille de seuils ce qui permet de calculer 
une pseudo courbe de pression capillaire. Cette approche a ete testee sur plusieurs 
6chanti lions de courbe de Pc connue, Les resultats montrent que Tobtention d*un bon 
accord avec les courbes de reference passe par une phase de calibration rigoureuse k 
realiser au cas par cas suivant la nature des roches etudiSes. Cette phase de calibration est 
30 n&essaire du fait de rincertitude au niveau : 



i er depot 
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© de la conversion signal RMN - distribution de taille de pore qui depend de la 
valeur de la relaxivitd de surface variable suivant les roches ; et 

* de la conversion distribution de tailles de pore - distribution de tailles de seuil 
qui depend de la nature de la roche et du processus de diagenese. 

-5 Cette approche_nIggt done pas recommandee dans un contexte predictif 

d'exploration. Elle ne serait de toute fagon pas applicable sur debris de forage. 

L* analyse d'image a aussi fait l'objet de travaux visant & l'obtention d'une courbe 
de Pc. Le milieu poreux est prealablement prepare sous forme d'une lame mince 
photographiee par microscopie electronique a balayage dite SEM pour "Scanning Electron 

10 Microscopy". L'image obtenue est alors analysee de maniere a determiner des parametres 
representatifs de la proportion et de la forme des vides par rapport k la roche. En 
particulier, il est possible de determiner une distribution de taille de seuil pour reconstruire 
une pseudo courbe de pression capillaire Pc. La principale limitation de cette methode est 
la nature bidimensionnelle (2D) de la lame mince tandis que la pression capillaire est par 

15 definition une propriete tridimensionnelle (3D). D'autre part, cette technique requiert un 
conditionnement assez lourd peu compatible avec Pobtention d'un resultat en leger diff6re. 
L' analyse d'image serait eventuellement applicable sur debris de roche mais n6cessiterait 
une calibration soigneuse pour acquerir une bonne predictibilite. 



20 chantier a des d6bris de forage mais e'est dans le but d'extraire le plus possible de fluide de 
forage de la roche pour minimiser les rejets polluants dans Fenvironnement mais aussi 
pour limiter les coQts en recyclant le fluide de forage r£cup6re. A notre connaissance, 
aucune mesure de centrifugation sur debris de forage n'a ete envisagee dans le but de 
determiner des proprietes capillaires. 

25 La methode selon Tinvention 

La methode selon Tinvention permet de determiner la courbe de pression capillaire 
de roches d'un gisement souterrain ipartir de mesures sur des debris ou fragments de 
roche qui y sont preleves (tels que des d6bris de forage), sur tout l'intervalle de saturation 
de ces roches, ceci dans un d61ai r6duit et & moindre cofit & partir de ces mesures. Elle 
30 comporte: 



A noter pour terminer que la technique de centrifugation est parfois appliquee sur 
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- une mesure de la permeabilite k des debris de roche ; 

- une mesure de la courbe de pression capillaire Pc en fonction de la saturation des 
debris de roche initialement satures d'un fluide en les soumettant a une centrifugation ; 
et 

- le param6trage dhone courbe de pression capillaire Pc satisfaisant k des relations 
empiriques dependant de parametres ajustables, que Ton contraint k s f ajuster k une 
partie asymptotique de la courbe capillaire mesuree par centrifugation, et k la valeur de 
la pernieabilite k mesuree sur les debris de forage, de fa§on k obtenir l'ensemble de la 
courbe de pression capillaire. 

Le parametrage de la courbe est effectue avantageusement en selectionnant par 
defaut un jeu des dits parametres permettant un calage sur la partie asymptotique de la 
pression. capillaire Pc aux faibles saturations, et par une modification de proche en proche 
des parametres de maniere que 1' estimation de la perm6abilit6 donn6e par une des relations 
empiriques utilisees, s'ajuste au mieux avec les mesures de permeabilite k realisees,, sur 
debris de roche et avec cette partie asymptotique. 

La mesure la permeabilite k des debris de roche est effectu6e par exemple k partir 
de mesures des variations de pression a l'interieur d'un recipient rempli d'un fluide 
contenant les d6bris de roche apr&s sa mise en communication durant un temps defini avec 
un reservoir de ce meme fluide sous pression, et a partir du volume effectivement absorbe 
par les debris de forage, et d'une modelisation de Involution de la pression ou du volume 
dans le r6cipient, a partir de valeurs initiales choisies pour les parametres physiques des 
d6bris de roche, que Ton ajuste iterativement pour que revolution mod61is6e de la pression 
s'ajuste au mieux avec revolution mesurSe des parametres physiques des debris de roche. 

La methode est avantageuse notamment en ce qu'elle fournit la pression capillaire 
des roches sur la base de simples debris de forage plus facilement disponibles et moins 
coflteux a obtenir .Lies resultats sont 6galement obtenus bien plus rapidement qu'avec des 
carottes. 
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1. Presentation sommaire des figures 

Les caracteristiques et avantages de la methode selon Tinvention, apparaitront plus 
clairement a la lecture de Ja description ci-apres d'un exemple non limitatif de mise en 
oeuvre, en se referant aux dessins annexes oil : 

- la fig.l montre revolution du signal RMN au cours de differents paliers de 



centrifugation fealis£s sur des debris de forage ; 



- la fig.2a montre les courbes de Pc obtenues sur une meme roche (B7) k partir de debris 
(cuttings) et d'une carotte ; 

- la fig. 2b montre la comparaison entre la courbe de Pc de reference mesufee sur carotte 
10 avec la courbe de Pc reconstruite sur tout 1'intervalle de saturation a partir des mesures 

de centrifugation et de la mesure ide la perm6abilit6 ; 

- les fig.3a et 3b montrent des resultats comparables k ceux montres aux Fig.2a et Fig.2b 
respectivement pour une autre roche GDV1 ; 

- les fig.4a et 4b montrent des resultats comparables a ceux montfes aux Fig.2a et Fig.2b 
15 respectivement pour une autre roche Rotl ; 

- les fig.Sa et 5b montrent des resultats comparables a ceux montres aux Fig.2a et Fig.2b 
respectivement pour une autre roche St Max ; 

- la fig.6 montre la comparaison entre les permeabilites de reference et les permeabilites 
mesurees sur debris de forage dans le cadre de la methode decrite dans deposee dans la 

20 demande de brevet FR 02/0023 precitee ; et 

- la fig.7 montre schematiquement rorganigramme de mise en oeuvre de la methode. 

Description detaillee de la methode 

La methode devaluation rapide d'une courbe de pression capillaire Pc a partir de 
deblais ou fragments de roche selon Tinvention est illustfee par la fig.7. La methode repose 
25 sur deux phases de mesure exp&imentales suivies d'une phase de param^trisation par 
reference avec des courbes connues, pour reconstruire la courbe de Pc sur tout 1'intervalle 
de saturation. Durant Tune des phases de mesure, on acquiert des donnees de pression 



capillaire Pc en centrifugeant les fragments de roche prealablement satures en eau. L'autre 
phase experimentale permet de calculer la valeur de permeabilite k de la roche a partir de 
la methode decrite dans la demande de brevet PR 02/0023 deja citee. La reconstruction de 
la courbe Pc sur tout 1'intervalle de saturation est effectuee en utilisant une forme 
parametree. Les parametres de la courbe sont determines de facon que la courbe s'accorde 
aux donnees de pression capillaire Pc obtenues experimentalement et avec la valeur 
estimee de permeabilite a partir de la courbe de Pc obtenue selon une m6thode connue telle 
que celle de Thomeer ou celle de Swanson Kamath qui seront rappelees plus loin, avec la 
valeur mesuree sur les fragments de roche. Plusieurs cas d' application sont proposes qui 
montrent le ties bon accord obtenu avec des courbes de r6ference sans proc6dure 
particuliere de calibration prealable. 

I) Mesure de Pc par centrifugation a partir des fragments de roche 

Pour la mise en ceuvre de la methode, il peut uuliser des moyens standards de 
centrifugation ou des moyens plus sophistiques avec un suivi automatique des volumes de 
fluide produits comme ceux decrits par exemple dans les brevets EP 603040 
(US 5,463,894), FR 2 763 690, PR 2 772 477 (US 6,185,985) ou PR 2 798 734 du 
demandeur. 

Les fragments de roche remontes au cours du forage sont prealablement nettoyes. 
avec des solvants dans un appareil de type Soxhlet, puis seches et satures en saumure 
30 g/1. Les fragments de roche sont ensuite 6gouttes dans un tissu humide de maniere a 
eliminer I'eau piegee entre les differents fragments de roche d puis introduits dans une 
cellule ou godet fixe a I'extremite d'un bras tournant. L'eau expulsee par centrifugation hors 
des debris de roche traverse une grille et est recueillie a la base du godet. L' acquisition des 
donnees experimentales se fait de la mSme maniere que dans le cadre de la centrifugation 
sur une carotte. Pour un palier de centrifugation (vitesse de rotation donnee), on mesure 
1'evolution de la production en eau jusqu'a ce que Ton n'observe plus de variations 
significatives, puis on augmente la vitesse de rotation pour commencer un nouveau palier. 

Comme le montre la fig.l on observe une diminution progressive du signal RMN et 
un decalage vers les faibles temps de relaxation T2 qui traduisent une desaturation du 
milieu poreux avec l'augmentation de la vitesse de rotation. Ce type de mesure prouve que 
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ron a bien un contact capillaire entre les debris de roche ce qui permet de mesurer la Pc 
par centrifugation a partir de debris de forage. 

Le volume d'eau produit au cours de Inexperience est converti en donnees de 
saturation a partir du volume d'eau initialement contenu dans les fragments de roche. Ce 
dernier est determine par pesee (diff6rence de poids des fragments de roche avant et apres 
farsatnr atiou) ou-directemeut-par-mesure-RMN. 



Les figures 2a a 5a montrent le resultat d'experiences realisees a partir de fragments 
de roche modeles, de taille 1 a 2 mm, fabriques au laboratoire a partir de roches de 
proprietes connues pour lesquelles on dispose d'une courbe de Pc mesuree de mamere 
classique par centrifugation d'une carotte. On. pent verifier que Von obtient une bonne 
correspondance avec la courbe de reference au niveau de la partie asymptotique (famle 
saturation en fluide mouillant). En revanche, on observe une difference importante pour les 
saturations en fluide mouillant plus elevees. On obtient done des resultats equivalents a 
ceux obtenus dans le cadre de mesures de porosimetrie utilisant le mercure sur fragments 
de roche et cela sans risque de pollution. 

La courbe de pression capillaire Pc mesuree n'est cependant representative que sur 
la partie asymptotique. Une procedure de reconstruction est done necessaire pour 6valuer 
le comportement de la courbe sur tout l'intervalle de saturation. 

II) Mesure de la permeabilite K des fragments de roche 

Pour mesurer la permeabilite des fragments de roche, on applique la methode 
decrite dans la demande de brevet precitee 02/02242. A cet effet, on plonge les fragments 
de roche dans une enceinte de confinement contenant un fluide visqueux. On met alors 
Venceinte en communication avec un recipient contenant ce mSme fluide sous pression, de 
maniere a comprimer le gaz piege dans les pores de la roche. Suivant un premier mode de 
mise en oeuvre. cette mise en communication pent etre tres breve et suivie apres un temps 
d'attente, d'une mesure de revolution de la pression dans renceinte. Suivant un aulxe mode 
de mise en ceuvre, la mise en communication est assez longue pour que Ton pmsse 
observer et mesurer la variation du volume effectivement absorbe par les fragments de 
roche. 



Ensuite, on mod£lise revolution de la pression ou du volume dans l'enceinte, k 
partir de valeurs initiales choisies pour les parametres physiques des fragments de roche, et 
on ajuste iterativement les valeurs des parametres physiques des fragments de roche de 
roche pour que devolution modelisee s'ajuste au mieux avec revolution mesuree du 
5 paramfctre physique dans 1 ' enceinte. 

Cette procedure donne d'excellents r6sultats. Les valeurs de permeabilite k des 
fragments de roche sont tout k fait conformes aux mesures de reference obtenues k partir 
* de carottes. 

HI) Reconstruction de la courbe totale de pression capillaire 

10 Dans cet te troisieme etape on fait la synthese des mesures precedentes ayant permis 

de construire la partie asymptotique de pression capillaire Pc et les mesures de la 
perm6abilite k, en se servant de relations empiriques, r6putees bien mod61iser; les 
param&tres physiques des roches. Par les publications suivantes : 

° Thomeer, J. H. M. : "Introduction of a pore geometrical factor defined by the capillary 
15 pressure curve", Trans AIME, vol , March, pp 73-77, 1960 ; et 

s Thomeer, J. H. M. : "Air permeability as a function of three pore network parameters", 
Trans AIME, vol , April, pp 809-814, 1983. J * 

sont connues des procedures devaluation de la permeability a partir d'une courbe de 
pression capillaire. La courbe de pression capillaire Pc est modelisee sous la forme 

20 P c = P d x exp( ^— ) ou : 

Ln ( Vb(Fc) ) 
V b (P./ 

o G est un parametre de forme pour prendre en compte la courbure de la courbe de 
pression capillaire (reliee a la forme de la distribution de taille de pore) ; 

° Pd> la pression de d£placement extrapolee k S Hg 6gale k zero ; et 

« le pourcentage de volume occupd par le mercure en fin d'exp&ience k pression 
25 capillaire infinie (egal k 4>xS Hg ), 



# 




les trois parametres 



10 

du modele etant relies a la perm6abilit6 par 1'expression suivante 



10 



k = 3.8068 xG -1 334 x 



P d J 



Dans la publication suivante : 
- Sanson. B.F. :» A simple cor relation between perm eability and mercury 
capillary pressures", JPT, vol , December, pp 2498-2504, 1981, 

n est P ropos6 egalement de cooler la valeur de la permeabilit6 avec la valeur 
maximale du rapport WP C ) sur la courbe de porosimetrie mercure. Ce point particulier 
correspond generalement au changement de regime qui se produit a la fin du regime de 
percolation, juste avant 1' augmentation significative de la pression capillaire. D'un point de 
vue capacite d'ecoulement de fluide dans le milieu poreux, ce point est particulierement 
important car il represente la tulle de pore pour laquelle Fensemble du reseau poreux est 
connects et qui controle done 1'ecoulement. La correlation la plus g£nerale foumie par 
l'auteur est donn6e par 1'expression (Swanson 1981). 



k =355x1 — ^ 



,a.ooa 



15 La publication suivante : 

- Kamath, J. : "Evaluation of accuracy of estimating air permeability from mercury 
injection data", SPE Formation evaluation, vol 7, 4, pp 304-310, 1992, 

porte 6galement sur revaluation comparative des approches par correlation empirique ou 
par modele physique pour determiner la valeur de la permeabilite a partir d'une courbe de 
20 porosim6trie mercure. Le meilleur accord est obtenu avec une nouvelle correlation basee 
sur la longueur caracteristique de Swanson 198 1 deja cite : 

i=413x4 8 ^ sik<lmD 
Jt=347x4 6 £ sik>lmD 
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Lmax etant d6fini de la maniere suivante : 



r«i»<s n 



t]>x>,x(100-S r ) 



et 



ioo-s r 



Avec : 

* o X : exposant traduisant la courbure de la courbe de pression capillaire (relive k la forme 
5 de la distribution de taille de pore). 

© P c : la pression de d6placement extrapolee a Sn g egale a zero. 

o S r : saturation r6siduelle occupee par le fluide mouillant (%). 

Pour parametrer la courbe de pression capillaire Pc empirique issue des^approches 
de Thomeer, Swanson ou Kamath, on la contraint k s'ajuster k la partie asymptotique 

10 obtenue par centrifugation durant la premiere etape de la m6thode Pour contraindre 
r ensemble de la coui'be, on se sert en outre de la valeur de la permeabilite k mesuree sur 
debris de forage durant la deuxieme etape de la methode, que Ton compare avec ce que 
donnent les relations empiriques. On modifie alors les parametres de la pression capillaire 
Pc jusqu'& ce que Ton honore k la fois le comportement asymptotique mesur6 et 

15 Festimation de la permeabilite ce qui permet de contraindre la courbe de pression capillaire 
sur tout Tintervalle de saturation Sw. 

Le processus d'inversion demarre avec un jeu de parametres par defaut qui 
permettent de caler le comportement asymptotique de la pression capillaire Pc aux faibles 
saturations en eau. On modifie alors de proche en proche ces param&tres (principalement la 
20 pression d* entree, not6 Pe ou Pd, et le facteur de forme, not^ X ou G) de maniere que 
Festimation de la perm6abilite donnee par une des relations precedentes soit en bon accord 
avec la mesure de perm6abilit6 realis6e sur debris de forage tout en conservant un bon 
accord avec les mesures de Pc aux faibles saturations en eau. 



1.1. Resultats obtenus 



25 Les figures 2b, 3b et 4b et 5b montrent la comparaison entre la courbe de Pc 

reconstitute suivant la proc6dure prec6dente et la courbe de Pc de reference obtenue sur 
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one carotte. On constate que, quel que soit People consider*, la methode de 
xecons.ru coon permet d'obtenir une evolution pertinente de la Pc sur tout rintervalle de 
saturation et en particulier aux fortes saturations en eau, quelle que soil la perm6ab,ht« de 
la roche. 





Nom 


K 

(md) 


G 


Vb 
(fraction) 


Pe 

(bar) 1 






Rotl 


150 


0.28 


0.89 


0.38 




GDV1 


195 


0.28 


1.00 


0.41 


B7 


780 


0.25 


0.92 


0.2 


StMax 


j 2000 


0.34 


0.92 


0.26 
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REVENDICATIONS 

1) Methode de detennination de la courbe de pression capillaire de roches d'un 
gisement souteirain a partir de mesures sur des debris de roche qui y sont preleves, 
comprenant une mesure de la permeabilite K des debris de roche, caract&risee en ce quelle 

5 comporte : 

- une mesure de la courbe de pression capillaire Pc en fonction de la saturation en 
soumettant a une centrifugation des d6bris de roche initialernent satur6s d'un fluide ; et 

- le parametrage d'une courbe de pression capillaire Pc satisfaisant a des relations 
empiriques dependant de param&tres ajustables, que Ton contraint & s'ajuster a une 

10 partie asymptotique de la courbe capillaire mesuree par centrifugation, et a la valeur de 

la permeabilite k mesuree sur les debris de forage, de fagon k obtenir l'ensemble de la 
courbe de pression capillaire. 

2) Methode selon la revendication 1, caracterisee en ce que le parametrage 
comporte la selection par defaut d'un jeu des dits parametres permettant un calage sur la 

15 partie asymptotique de la pression capillaire Pc aux faibles saturations, et une modification 
de proche en proche des paramfetres de mani&re que r estimation de la permeability donn6e 
par une des relations empiriques .utilisees s ! ajuste au mieux avec les mesures de 
permeabilite realisees sur debris de roche et avec la dite partie asymptotique. 

3) M6thode selon la revendication 1 ou 2, caract6ris6e en ce que Ton mesure la 
20 permeabilite k des d6bris de roche k partir de mesures des variations de pression h. 

I'int6rieur d'un recipient rempli d'un fluide contenant les d6bris de roche apres sa mise en 
communication durant un temps defini avec un reservoir de ce meme fluide sous pression 
et k partir du volume effectivement absorb^ par les d£bris de forage, et d'une modelisation 
de Involution de la pression ou du volume dans le recipient, k partir de valeurs initiales 
25 choisies pour les parametres physiques des d6bris de roche, que Ton ajuste iterativement 
pour que Involution modelisee de la pression s'ajuste au mieux avec revolution mesuree 
des paramfetres physiques des debris de roche. 
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